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SCHEMES FOR MULTI-DIMENSIONAL ISAACS EQUATIONS. 
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Abstract. Recently, Krylov, Barles, and Jakobsen developed the theory for 
estimating errors of monotone approximation schemes for the Bellman equa- 
tion (a convex Isaacs equation). In this paper we consider an extension of 
this theory to a class of non-convex multidimensional Isaacs equations. This 
is the first result of this kind for non-convex multidimensional fully non-linear 
problems. 

To get the error bound, a key intermediate step is to introduce a penal- 
ization approximation. We conclude by (i) providing new error bounds for 
penalization approximations extending earlier results by e.g. Benssousan and 
Lions, and (ii) obtaining error bounds for approximation schemes for the penal- 
ization equation using very precise a priori bounds and a slight generalization 
of the recent theory of Krylov, Barles, and Jakobsen. 



1. Introduction 

In this paper we will study error bounds for approximation schemes for a class 
of non-convex multidimensional Isaacs equations. To be precise, we will consider 
the following (one-sided) obstacle problem 

(1.1) mm{F{x,u,Du,D'^u),u- g} ^0 in M^, 
where g is the obstacle (u > g), and F is given by 

(1.2) F{x,r,p,X) = snp {-ti[a°'{x)X] + b°'{x)p + r{x,r)} . 

Here A is a compact metric space, a is a positive semidefinite matrix, / is strictly 
increasing in r, and the data is a least bounded and uniformly continuous. Precise 
assumptions will be specified later. This equation is non-convex because of the 
min / sup form, but also because the / term may be non-convex in r. It may also 
be degenerate since a may vanish at certain x and a. Under the assumptions we 
will use, this equations can always be rewritten as an Isaacs equation, 

inf sup l-tr \a"'-'^(x)D^u\ + h°''^(x)Du + c"''3(x)w + n^ix)] = 

in R-'^, for suitably defined a, 6, c, and /. It is well-known that such equations do 
not in general have smooth solutions. 

The above problem Hl.l|l is also called a variational inequality, and such prob- 
lems occur in many applications and have been studied intensively for a long time. 
The classical theory for variational inequalities (see e.g. j^llHI) studies weak or 
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variational solutions and uses cither PDE techniques and Sobolev space theory or 
probabilistic techniques using the connection with optimal stopping time problems. 
In this paper we will (mostly) consider viscosity solutions which is a weaker and 
more recent concept of solutions. We refer to (Sj (lO, for a general overview of the 
viscosity solution theory, and to e.g. [JSJ for analysis and applications of 

obstacle problems in the viscosity solutions setting. We mention in particular the 
many applications in finance, like e.g. the pricing problem for American options 

m- 

In the viscosity solution setting the first results on error bounds for monotone 
schemes were obtained by Crandall and Lions for first order equations. This case 
has later been studied by many authors. Only recently did Krylov OH] obtain 
the first results for second order fully non-linear equations (the convex Bellman- 
equations), and these results were then extended and improved by Barles and Jakob- 
sen O 0] . We refer to the recent paper ^ for the best results available at the 
present time. All these results concern the convex Bellman equation. In the non- 
convex fully non-linear case, there are to the best of the author's knowledge no 
results in the multi-dimensional case. The only non-convex result we know about 
applies to one dimensional problems T4'. 

In this paper we will give error bounds for general monotone approximation 
schemes for the non-convex multi-dimensional problem Hl.l|) . We will use the fol- 
lowing abstract notation for such schemes. 



(1.3) min {S{h, x, Uh{x), [uh]x)] Uh[x) - g[x)} = 



7,N 



where S is loosely speaking a consistent, monotone, and uniformly continuous ap- 
proximation of F in The approximate solution is Uh, [uh]x is a function 
defined from Uh, and the approximation parameter is h. This notation was in- 
troduced by Barles and Souganidis j^l to display clearly the monotonicity of the 
scheme: S is non-decreasing in Uh and non- increasing in [w/ija;. Typical approxi- 
mations S that we have in mind are certain finite difference methods (FDMs) pi] 
and so-called control schemes (Zj. In Section [3 we will explain the notation for a 
concrete FDM. 

To get an idea of our results, we will now consider an explicit ID problem: 



mm sup 



|-a"A2w" + /"(a;,u)};M-5(a:)| =0 in R\ 
We approximate this problem using a monotone FDM, 

mini snp a" Aluh + f"ix,Uh)y,Uh- gix)i =0 in R\ 



^aeA 

where 

2 ^ (/)(X + h)- 2(j){x) + (^{x - h) 

^hH^) = • 

Under suitable assumptions on / and g our results (cf. Proposition 15.1(1 give the 
following error bound, 

\\u-uh\\l^ < Ch'/^ 
If we take slightly stronger assumptions on the obstacle g, we get 

\\u-Uh\\L^ <Ch'/\ 
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The same results hold for problems in arbitrary space dimensions, see Sectional In 
the convex case {g = —oo) the corresponding rate is 1/2 jl5| . 

In the (much) more difficult case when a also depend on x, the rate in the convex 
case is a least 1/5 In this paper we will not be able to handle FDMs when the 
coefficient a depends on x. Results in this direction can be obtained by combining 
the methods of this paper with those of . But since the arguments become much 
more involved in this case, we have chosen to omit it. However note that the results 
in this paper applies to the so-called control schemes in the general case (when a 
depends also on x), see Section |S1 For a better discussion of this point we refer to 

Ei 

Let us now try to explain how we get our error bounds. As a key intermediate 
step we introduce the following penalization problem, 

(1.4) F{x,Ve,Dve,D\>,)^^{v,^g)~ in R^, 

where (•)^ = — min(-,0). Under suitable assumptions on the data, it is possible to 
show that the solution of the penalized problem H1.4|l converges monotonically to 
the solution of the obstacle problem Hl.l(l as £ ^ 1_ . In this paper we prove new 
error bounds for this convergence using easy comparison arguments. 

The next step is then to consider the approximation scheme associated to ()1.4|l 
via ((T3|l . 

(1-5) S(h,x,Vh,^{x),[vh,e]x) = - {vh,6 - gV in M^. 

£ 

Again we prove that Vh,e converge to Uh the solution of (|1.3|l with a given error 
bound. This argument is completely similar to the one mentioned above in connec- 
tion with H1.4I) . 

The third and more difficult step is to obtain error bounds for convergence the 
solution Vh,e of (|1.5|l to the solution of (|1.4|l . To get this result we use a slight 
extension of the arguments in O I15| . What is new here is that the equation need 
not be convex in the zero-order term (the u term), as is the case for (|1.4ll . 

As a finial step we combine the previous steps to get the full error bound via the 
triangle inequality, 

\\u - < ll" - V^Wl-- + \\Ve - Vh.dlL'^ + \\vh,e " UftlU^- 

The right hand side will depend on h and £, and the result follows after a mini- 
mization over £ > 0. Warning! This last step is only possible to perform if the 
bound on \\ve — Vh,e\\L°° does not depend on £ in a too singular manner. Note that 
some coefficients in H1.4|l and ()1.5|l depend on l/£, and that naive computations 
would lead to a priori bounds on the solutions that also depend on l/£. With such 
bounds, we would not be able to prove any error bounds. For our purpose, we need 
and prove more precise a priori bounds than can be found in the literature. 

Let us now return briefly to the penalization problem (|1.4|l . Usually it is easier to 
obtain existence of solutions of the penalized problem and than of the corresponding 
obstacle problem. The limit procedure (£ 0) then gives existence of also for the 
obstacle problem. We refer to Bensoussan and Lions [S] for the classical theory 
and to Amadori pP for a viscosity solutions approach. Error bounds exist in the 
classical case. E.g. in [HJ p. 197] the following bound is proved, 

(1.6) \\v,-u\\w^.2 <Ce^^^, 
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in the case when F in is hncar, uniformly elliptic, and in divergence form. In 
this paper we prove under suitable assumptions that 

(1.7) \\v,-uh^<Ce'^\ 

and under sightly stronger assumptions on g we get Hfg — uIIloo < Ce. These results 
applies to very general equations, see Section |21 to all kinds of weak solutions as 
long as the comparison principle holds, and even to monotone schemes like (|1.3f) 
and (|1.5I) . To the best of the author's knowledge this result is new, even in the 
linear uniformly elliptic case e.g. under the assumptions leading to H1.6|l . Also note 
that l|1.7|l does not follow from (|1.6(l except in one space dimension (by Sobolev 
embedding) . 

Let us now introduce some notation: We will use the following (semi) norms, 

|/|o = esssup^gK«|/(a;)|, [/]^ = esssup^,j^gRiv^^^^^— |/U = l/|o + [/l^, 

where / : M*^ is a function and /i G (0,1]. The same notation will be 

used for vector and matrix valued functions /, in which case |/| is interpreted 
as a vector and matrix norm respectively. L°°(E^), C(M^), Ch(M^), C°''"(M^), 
At G (0, 1], C'=(M^), ken, denote the spaces of functions f : ^ R that are 
bounded, continuous, bounded and continuous, have finite norm \f\p_, and are k- 
times continuous differentiable respectively. Furthermore, W^'^, Wq'^, W'^''^, and 
W^'°° = C*^'^ are standard Sobolev spaces. The space of real symmetric N x N 
matrices are denoted by S^, and X > Y in E>^ will mean that X — Y is positive 
semi-definite. Finally, by D'^cj) we mean the vector of fc-order partial derivatives of 
a function (f>. 

The outline of the rest of this paper is as follows: In the next section, we treat 
the penalization problem (|1.4|l . We state and prove a very general error bound and 
compare it with classical results by Benssousan and Lions. In Section |3| we obtain 
error bounds for equations that are non-convex in the 0-th order term. These results 
are of auxiliary nature and are needed in Section 01 In this section we state and 
prove our main result, an error bound for ^l.l^ . Then in Sectional we apply our 
main result to obtain error bounds for a FDM and a control scheme. Finally there 
is an Appendix containing some technical a priori estimates. 

2. The penalization method 

In this section we will use comparison arguments to derive new error bounds for 
the convergence of the solution of the penalization problem (ll.4|l to the solution of 
the obstacle problem 1)1.1(1 . Here we will no longer assume (ll.2|l . in stead we will 
allow for very general structure of F: 

(CI) (Comparison) The equations 1(1. 1() and (|1.4() satisfy the comparison principle 
for the class of weak solutions under consideration. 

(C2) (Monotonicity) Let X, F e §^,p, x € R^, r, s e K. If X > F and r < s, then 

F{x,r,p,X) <Fix,s,p,Y). 

(C3) (Regularity) One of the following statements hold: 
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(i) g e CO'i(M^), \D'^g-\ < C, and for every x e and (j) e C^{R'^) 
satisfying \(j)\o,i + \D^(j)-\o < R, 

F{x,4>{x),D4>{x),D''cp{x)) < Cr. 

(ii) g e C"^i(R^), and for every x G and (j) G C^{R^) satisfying \(t>\o,i < R, 

F{x, <t){x),D<t>{x),D''cp{x)) < Cr{1 + \D^(j)-\o)- 

(iii) g e C°''^(R^) for some £ (0, 1), and for every a; € R^ and G C'^{R^) 
satisfying |(/)|o < R, 

F{x, (f^ix), D(f>{x), D^cPix)) < Cnil + \D^\o + \D^- |o). 

Assumption (CI) is not very precise. In applications we need to specify both 
the notion of weak solutions and "boundary conditions" at infinity. Assumption 
(C2) says that F is "proper" in the terminology of the User's Guide 'S', and implies 
that F is degenerate elliptic. Assumption (C3) gives regularity assumption on the 
obstacle g and corresponding (local) boundedness assumptions on F . Assumption 
(C3) can be generalized to allow for super- linear growth in \X^\ and \p\. This 
would affect the rates obtained and will not be considered here. 

These assumptions are satisfied by a very wide class equations and with different 
concepts of weak solutions. In the viscosity solutions setting (the weakest notion 
allowed here), we will just mention that the above assumption hold for the Bellman 
equations from stochastic control ^01 and the Isaacs equations from stochastic 
differential games ^I] under natural assumptions on the data. We refer to the 
User's Guide [H] for many more viscosity solution examples. Typical "boundary 
conditions" would be to assume bounded solutions or linear growth at infinity. We 
can also consider variational solutions [H] whenever it makes sense to do so. In 
this case all point-wise inequalities have to be interpreted in the almost everywhere 
sense. 

The main result in this section gives both the convergence and the rate of con- 
vergence for the penalization problem. 

Theorem 2.1. Assume (CI) - (C3) hold and u and are solutions of and 
H1.4I) (we do not assume p.2|l .^). Then if \D'^g~\Q < oo (case (i)) 

Q<u-Ve<Ce in R^, 

otherwise (g G C^'^^{K.^) - cases (ii) and (iii)) 

< u - We < Ce'^/^ in R^, 

where the constants C only depend on g and Cr from (C3). 

Remark 2.1. The rates depends only on the regularity of the obstacle, and not on 
the regularity of the solution. Even if the solution u is only Holder continuous, we 
still get rate 1 if |-D^.g^|o < oo. For many other types of approximation schemes 
the rates depends directly on the regularity of the solution, see e.g. |3| (FDMs) and 
jl7| (vanishing viscosity method). 

To the best of the our knowledge, this is the first time the penalization error has 
been estimated for degenerate equations, and the above result seem to be new even 
in the linear uniformly elliptic case (see below). 

Before giving the proof of Theorem l2.1l let us briefly consider the linear uniformly 
elliptic case. Here L^/H^'^'^-estimates on the penalization error are classical 0. We 
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will state a typical such result, so that the reader can compare it with the one we 
have obtained. In our notation: 

(2.1) mm {-Au + f{x);u - g{x)} = in Q, 

(2.2) -Av, + fix) = ^{v,-gix))- in n, 

where f2 is a smooth bounded domain and A is a linear elliptic operator in divergence 
form, 

A(j){x) := dxi{aij{x)dxj(l>) + bi{x)di(f)- X4>. 

The summation convention is used, A > 0, and ellipticity means ^iaij{x)^j > 
for some a > and every ^ S M.^ . The concept of solutions is that of variational 
(weak) solutions belonging to Wg^'^(ri), see |S] for the exact definitions. Typical 
assumptions on the data are 

(D) a^jeL^{n),heW^'°°{n)J £L^{n),geW^''^{n), AgeL^n), 

where a = {aij)ij,b — (6i)i, and the error bound obtained is the following 6, p. 
197]: 

Proposition 2.2 (Bensoussan & Lions). Assume (D) holds, A > large enough, 
g\dn > 0, and u and solve (|2.1|l and (|2.2|l . Then 

||w-^'e||wi.2(o) < Ce^/^. 

Note that in this theorem we need control over the second derivatives of the 
obstacle g [Ag g essentially means that g E W'^''^), while in our result we only 
need to control the first derivative of g (say g G W^'°° = C^'^). Furthermore, we 
may use Theorem 12 .11 to get a new error bound in this case. Comparison principles 
for H2.1|l and H2.2|l are essentially given by Theorems 1.2 and 1.4 p. 192 and p. 198 
in [5], so if (? G W^'°°{^1) we can conclude by Theorem 12 . II that 

\\u-Ve\\L^^n)<Ce'^\ 

Here u and are variational solutions of H2.1|l and H2.2|) . This result does not 
follow from Proposition 12 . 21 unless is a domain in M}. 

The proof of Theorem 12.11 We give a series of simple lemmas that leads the 
way to the proof of Theorem l2.1l We start by a preliminary error estimate: 

Lemma 2.3. Assume (CI) and (C2). Let u and solve Hl.l|l and ()1.4|l . Then 

Q<u-v^<\{v^-g)~\^ in M^. 

Proof. First we check that by monotonicity in r (C2), 

Ve + \ {Ve - 9y\o 

is a supersolution of 1)1.1(1 . The comparison principle for 1(1. l|l then yields the second 
inequality. Similarly, the first inequality follows since is subsolution of 1(1.1(1 . □ 

Now we will estimate \{v^ — 5)^|o: 

Lemma 2.4. Assume (CI) - (C3) hold and g e C^(E^). According to (C3) define: 
Case (i): K := Cr with R = |g|o + \Dg\Q + \D'^g-\Q. 
Case (a): K := Cr{1 + |-D^ff"|o) with R = |g|o + \Dg\o. 
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Case (ill): K Cr{\ + \Dg\Q + \D'^g-\o) with R = |g|o. 
Let Ve he the solution of H1.4() . Then 

-£K<v,-g m M^. 

Proof. The result follows from the comparison principle since 

g-sK 

is a (classical) subsolution of (|1.4() . □ 

Since we did not assume that g is smooth, we need an approximation result. 
Let ps be the standard moUificr, ps{x) = 1/5^ p{x/5) where p is a smooth positive 
function with mass one and support in the unit ball. Let gs = g * ps and denote by 

and vl the solutions of Hl.l|l and H1.4|l when gs has replaced g: 

(2.3) TtiwL{F{x,u^,Du^,D'^u^)-u^ ~gs}^0 in M^, 

(2.4) F{x,vlDvlD\l) = -^{vl-gsr in M^. 
We have the following bounds on u — and — v^: 

Lemma 2.5. Assume (CI) and ( C2), and let u, , and be solutions of 
(lO) . (O, and (E3|l- Then 

\u-u'\o + \v,~v% < \g~g'\o. 

Proof. We only prove the w-result, the proof of the u-rcsult is similar. (If we knew 
a priori that ^ u, the w-result could be obtained by going to the limit in the 
w-result.) Let K := \g — g^\o and define 

w=^(x) =v^,{x)±K. 

The result follows by the comparison principle for l|1.4() since z«+ and w~ are super- 
and subsolutions of (|1.4|) respectively. 

Let us prove that is a supersolution of (|1.4|) . the subsolution part is similar. 
First observe that by (C2) 

F{x,w+,Dw+,D^w+) > F{x,vl,Dvl,D'^vl). 

Then observe that by the definition of K , 

-{w+ - g)- = -{vi + K- g)- > -{vl - gs)- . 

Since is a supersolution of (|2.4|) , the above observations show (at least formally) 
that 

F{x,w+,Dw+,D^w+) - i(w+ - .g)- 
> F{x, vl Dvl D^vl) - -^{vl - gs)- > 0. 

The proof is complete since all the above computations easily can be seen to hold 
in the weak/ viscosity sense. □ 

Now we can give the proof of Theorem 12. II 
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Proof of Theorem \2.1\ First we consider the solutions and of (|2.3(l and H2.4|l . 
By Lemmas O and O we have 

~vl< Ke in M^, 

where K is defined in Lemma 12.31 Since 

u~v,^{u~u^) + [u^ - vl) + {vi - V,), 

Lemma 12.51 and Holder continuity of g lead to 

u-v,< 2[g]^,5^' + Ke. 

If IffaT^-lo < oo then K is independent of 5 and we can send 5 — > 0, leading to 

u — Ve < Ke. 

Otherwise, by Holder continuity oi g, K = C5^~'^ , and minimization w.r.t. 5 yields 

u~v,< Ce^/^ 

The lower bound follow from Lemma 12.31 □ 

Remark 2.2. The procedure used in the above proof is very general, and works for 
any problem satisfying the assumptions corresponding to (CI) - (C3). E.g. one 
could consider boundary value problems where one would find that the estimates 
in Theorem 12 . II still holds. In the next section we will even see the method applied 
to an obstacle problem for an approximation scheme ('Lemma l4.1|l . 

We end this section by indicating an alternative approach. Remember that 
the lower bounds in Theorem 12.11 follow from the comparison principle since the 
solution Ve of 1)1.4(1 is a subsolution of p. 1(1 . To obtain the upper bounds, we need 
in some sense to show that is an approximate supersolution of 1(1.1(1 . Observe 
that formally 

m.iii{F[ve\]v, - g} > -e{F[v,])^. 

This follows since 

= F[v,] - -{v, - g)- = min \ F[v,]; F[v,] -^-{v,^ g) 
e ye 

implies that 

= Ta:\TL{F[v^];eF[v^] +v^- g) < Tah\{F[v^]\Ve - g} + e[F[ve]y . 

This is vanishing viscosity(!), and we should already guess that the error should be 
Ce^/^ when solutions are Lipschitz continuous. An easy way to get this result is 
the continuous dependence approach of ^1 El which leads 

u~v,< C"e(F[«,])+ " - C(l + \Dv,\o)e^/\ 

where u is the solution of 1(1. 1() . The loss of rate is caused by being only Lipschitz 
continuous while is a second order operator. On the other hand, if |Z3^u^|o < oo 
then wc would have the full rate: 

u^v,< = C{1 + \Dv,\o + \D\-\o)e. 

In Theorem 12.11 this last estimate is proved under the much weaker assumption 
\D'^g-\o < CO. 
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3. Monotone Approximation Schemes - Preliminaries. 

In this section we will give a slight generalization of the results of ^| EH El ■ 
We prove error bounds for monotone approximation schemes for equations that 
have possibly non-convex dependence on 0-order terms. These results will then be 
used in Section^ to obtain rates for the more difficult obstacle problem Hl.l|l and 

Consider the following equation: 

(3.1) F{x,u,Du,D^u) = in R^, 

where F is given by (|1.2(l in the introduction. We make the following assumptions: 

(Al) a" = ^a°'a°'^ for some N x P matrix a, and there is a C independent of a 
such that |(t"|i + < C. 

(A2) There are A, A > such that for every x G M^, a £ A, r, s E R satisfying 
r > s, 

A(r - s) < rix, r) - /"(x, s) < A(r - s). 

Furthermore, /"(-jO) is bounded uniformly in a, and for every x,y £ M^, r € R, 
and a G A 

\rix,r)~r{y,r)\<C{l + \r\)\x-y\. 

Remark 3.1. The first part of assumption (A2) implies that / is Lipschitz and 
strictly increasing in r. The second part implies that / is bounded and Lipschitz 
in x for fixed r. If (Al) and (A2) hold and e > is fixed, the penalization scheme 
(|1.4|l can be rewritten in the form H3.1|) by redefining /"(x, r) to be 

/"(x, r) H — min{r — g{x); 0} . 

This new function than satisfies (A2) with new constants A, A + i, and C. 

Existence, uniqueness, and regularity follow from standard viscosity solutions 
arguments. The results parallels the one mentioned in Section ^ and we state 
them without proofs: 

Lemma 3.1. Assume (Al) and (A 2) hold. Then there is a unique bounded Holder 
continuous viscosity solution u of (13.11) . Furthermore, if X is big enough (compared 
to [a]i and [b]i), then u is Lipschitz continuous. 

Using notation from the introduction, we may write an approximation scheme 
for l|3.1|l in the following way 

(3.2) S{Kx,Uh{x),[uh]^)^Q in R^. 
We require S to satisfy: 

(SI) (Monotonicity) For every /i > 0, x € R^, r € R, to > and bounded functions 
w, V such that u <v m R^, the following holds: 

S{h, x,r + ni,[u + m]x) > Am + S{h, x, r, [v]x), 

where A > is given by (A2). 
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(52) (Regularity) For every h > and cf) e Ch(M^), x ^ S{h,x,<j>{x),[<j>]x) is 
bounded and continuous in M.^ and the function r i-^ S{h,x,r, [<j>]x) is uniformly 
continuous for bounded r, uniformly in a; G M.^ . 

(53) (Consistency) There exists integers n, ki > 0, constants Ki > 0, i ~ 1,2, ... ,n 
such that for every smooth (f>, h > 0, and x e M^: 

n 

\F{x, ct^{x),Dct^{x),D^c^{x)) - S{h, x, 4>{x), [4>]x)\ < ^ K,\D'q^{x)\h''\ 

1=1 

Condition (SI) and (S2) imply a comparison result for bounded continuous so- 
lutions of 1221) (cf. 3 ): 

Lemma 3.2. Assume (SI), (S2), and u,v e Cfc(M^). If S[u] < and S[v] > in 
R^, then u < v in R^. 

We proceed with obtaining an upper bound on the error for the scheme H3.2|l . 
In order to do so we will consider the following auxiliary problem: 

(3.3) su-p F{x + e,u^{x),Du^{x),D^u\x)) = in R^, 

\e\<S 

where 6 > 0, and with u being the solution of (|3.1|l . 

F{x,r,p,X) := sup {-tT[a"{x)X] + b°'{x)p + Xr - Xu{x) + f°'{x,u{x))} . 

Actually this is a problem of the same type as 1)3. 1() so well-posedness follows in the 
same way. At this point we assume the following: 

(A3) Let u and denote the solutions of H3.1|l and H3.3|l . There is a constant 
A' > independent of S such that 

\u^\i + ^\u^u^\o<K. 



Remark 3.2. Assumption (A3) follows from assumptions (Al) and (A2) if A is big 
enough. After observing that (|3.1|l can be written as an Isaacs equation, this follows 
from Lemmas lA. II and lA.2l in the Appendix. In the case that A is not "big enough" 
things are a little bit more complicated, we refer to 0^ for this case. 

Now we are in a position to derive an upper bound on the error for the scheme 
dial- 
Theorem 3.3. Let (Al) - (A3), (SI) - (S3) hold, let u be the viscosity solution 
of (I3.1|l , and let Uh be a solution of the scheme (|3.2|l . Then if h > is sufficiently 
.small, 



u — Uh < Ch'^ in 



wherej:= min and C < ^(V"_, AT^ + 2(2A + A)). 

Proof. 1) We start by showing that us := ps * u^ is a subsolution of 
(3.4) F{x,w,Dw,D^w) ^0 in R^, 

where ps is the mollifier defined in Section 13 By (A3) 

F{x + e,u\x),Du\x),D^u\x)) <0 in R^ 
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for every |e| < S. Hence for every |e| < S, u^{x — e) is a subsolution of H3.4|l . 
Then us is also a subsolution of (|3.4(l since it can be viewed as the limit of convex 
combinations of subsolutions u^{x — e) of the convex equation (|3.4|) . we refer to the 
Appendix in jSj for the details. 

2) Us is an approximate subsolution to the scheme (|3.2|) . By properties of mol- 
lifiers and (A3), us is smooth and satisfies 

5'~^\D'us\o + [25)-^\u-us\o<K. 

So by (A2) and the definition of F, for every x e R^, 

F{x, u\x),Du\x), D^u\x)) < sup \X{us - it) - /"(a:, us) + /"(x, u)\ 

< 2K{\ + K)5. 

Consistency (S3) then leads to 

n 

S{h, y, us{y), [us]y) < K,5^''h!'' + 2K{\ + k)6 =: C. 

2 = 1 

3) By (SI), Us — C/X is a subsolution to the scheme (|3.2|) . By comparison. 
Lemma 13.21 we have 

us-Uh< C/X in M^. 

4) Combining the above estimates yields 

u ~ Uh = u — Us + Us ~ Uh < 2KS + C/X in . 
Now we can conclude by choosing 

S= maxl/i*^-/*}. 

i:Ki>0 

□ 

Remark 3.3. If we replace sup by inf in equations H3.1|l and (|3.3|l . a similar argument 
would lead to a lower bound of the error: —Ch"' < u — Uh- 

From this remark it is clear that we have the full result for semi-linear equations 
(see also [21 for the linear case): 

Corollary 3.4. Assume l|3.1|) is semi-linear, i.e. that A is a singleton. Let (Al) 
- (A3), (SI) - (S3) hold, let u he the viscosity solution of (|3.1|l . and let u^ be a 
solution of the scheme (|3.2(l . Then if h > is sufficiently small, 

\u - Uh\a < Ch'', 

where 7 and C are defined in Theorem \3.3[ 

Following the ideas in , we proceed to have obtain the full result in for 
more general situations. Let S denote the scheme S when it is applied to equation 
F[w] =0 where F is defined just after 1)3. 3|l . and consider 

(3.5) sup S{h,x + e,ui{x),[ui]^) in R^, 

|e|<<5 

and the assumption analogous to (A3): 
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(54) Assume Uh and uf^ are solutions of (|3.2(l and (j3.5|l . and there is a constant 
K' > independent of S such that 

\ui\, + ^\uh~u%<K'. 

In addition we need the foUowing assumptions of S: 

(55) (Convexity) For any v € C^'^{M.^), /i > 0, and a; £ 

/ S{h, X, v{x - e), [v{- - e)]x)ps{e)de > S{h, x, {v * ps){x), [v * ps]x)- 

(56) (Commutation with translations) For any h > small enough, < (5 < 1, 
y e R^, teR,ve C'biR'^) and |e| < S, we have 

^(/i,y,t, H^_J = 5(/i,y,t, [«(--e)]J). 

Remark 3.4. While (S5) and (S6) are not very restrictive, (S4) is. This assumption 
is satisfied for control schemes in general and for FDMs when the coefficients 
multiplying second order derivatives are constants (Section EJ. Note that (S4) is 
not assumed in Corollary 13. 41 

It is clear that by repeating the arguments in the proof of Theorem 13.31 with 
the schemes (|3.2ll and (|3.5|l taking the role of the equations (|3.1I) and H3.3() . we 
obtain a lower bound on the error —CK^ < Uh — u. We refer to 3 for more details. 
Combining this result with Theorem 13 . 31 then yields the main result in this section. 

Theorem 3.5. Let (Al) - (A3), (SI) - (S6) hold, let u he the viscosity solution 
of (|3.1|l , and let Uh be a solution of the scheme l|3.2|l . Then if h > is sufficiently 
small, 

\u - Uh\a < Ch'^ , 

where 7 is defined in TheoremXT^ and C < ^» + 2(2A + A)). 

The results in this section generalize slightly the results in |19[ I2(JI O which 
consider pure convex or concave equations. Here we allow non-convexity (non- 
concavity) in the 0-th order terms. 

4. Monotone Approximation Schemes - The Main Result. 

In this section we will see how to use the results of the previous two sections to 
obtain error bounds for monotone schemes H1.3|l for the non-convex problem 1)1.1(1 . 
In Section[51we give examples of such schemes. We assume that S in (|1.3|) satisfies 
assumptions (SI) - (S3) of Section |31 

First we consider the penalization problem corresponding to 1(1. 3|l . namely prob- 
lem ()1.5|l in the introduction. This scheme is also an approximation scheme for the 
penalization problem ((1.4|l . Note that 1(1.3(1 and ((1.5(1 themselves satisfy assump- 
tions (SI) and (S2) (when S is appropriately redefined) and hence the comparison 
principle. Lemma [3. 21 holds also for these schemes. 

We start by obtaining the rate of convergence for vf^ — ^ Uh- It is not difhcult 
to see that this result is a consequence of the procedure given in Section |21 if we 
can prove that assumptions corresponding to (CI) - (C3) hold for S. Because (SI) 
and (S2) imply comparison, they already imply assumptions corresponding to (CI) 
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and (C2). But it turns out that (S3) is not sufficient to have the assumption corre- 
sponding to (C3) because it involves derivatives of higher order than two. We need 
to assume "(C3)": 

(S7) (Regularity) One of the following statements hold: 

(i) g G C"3'i(R^), \D^g-\o < C, and for every x e M.^ and (f> E C^{M.^) 
satisfying \<j)\„^i + \D'^<j)-\„ < R, 

S{h,x,(l){x),(f>) < Cr. 

(ii) g e CO'i(K^), and for every x € and </> G C^{R^) satisfying \(f>\o,i < R, 

S{h,x,cl){x),^) <Cr{1 + \D^(I)-\o). 

(iii) g e C"'^(R^) for some fi G (0, 1), and for every a; € and (/) G C^{R^) 
satisfying \(j)\o < R, 

S{h, X, 0(x), </.) <Cr{1 + \D<j>\o + \D^-\„). 

Remark 4.1. This assumption holds for most reasonable schemes (|1.3ll when (S3) 
also holds, e.g. for the finite difference method (|5.1() below. 

By the method of Section 12 we have the following result: 

Lemma 4.1. Assume (SI), (S2), (S7) hold and and v^^e a'^e solutions of p. 3(1 
and 1)1. 5f) . Then if\D^g^\Q < oo (case (i)) 

0<Uh- Vh.e < Ce m R^, 

Otherwise (g € C'^'^(R"'^) ~ cases (ii) and (Hi)) 

0<Uh- Vh,e < Ce^'l'^ m R^, 

where the constants C only depend on g and Cr from (C3). 

Now to obtain results for the scheme p.3|l . we may use the following diagram: 

min{F[u]; u — <?} = ^ '■ min{S'[u;i]; u — 5} = 

0<M-i>E<Ci(e) 

F[Ve\ = \{Ve - gy < IhcoremLLa S[vh,e] = jivh,e - gT 

The main result of this paper is the following: 

Theorem 4.2. Let (Al), (A2), (SI) - (S7) hold with A > aup^{[a"]l + [b"]i) m 
(A2) and K' independent of e in (S4), let u he the viscosity solution of 1(1. l|l with 
F defined in ((1.2(1 . and let he a solution of the scheme p. 3(1 . Then if h > is 
sufficiently small, 

\u-Uh\o<Ch-'^\ 
If in addition \D'^g^\o < 00, then 

\u~Uh\o<Ch-'^^. 

Here 7 is defined in Theorem \S.!A and the constants C are independent of h. 



Theorem im 



0<Uh~Vh,s<Ci(e) 
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The assumption on A may be relaxed to simply requiring A > 0. This will 
influence the rates and complicate the arguments, see |H] for a discussion. See also 
Remark 13.21 

Outline of proof . 1) By Lemmas lA.ll and IA.2I in the Appendix (see Remark l3.2|l . 
assumption (A3) is satisfied for H1.4(l with K independent of e! Note that K' is 
assumed independent of e. 

2) By Theorem ESI with A replaced by A + i, 

\ve-Vh,e\o<C{l + ^)h^ in M^, 

where v^.e solves (|1.5f) and C is a constant independent of e. 

3) The result now follows from the triangle inequality, part 2) , Theorem 12.11 
Lemma l4.ll and a minimization in e (see the above diagram). □ 

If F is concave instead of convex so that the obstacle problem (|1.1() is concave, 
then we obtain better rates using directly Theorem 13. 51 

This was essentially the case considered by [3 . Theorem 14.21 is the first result 
for multi-dimensional non-concave/non-convex equation. 

5. Applications 

5.1. A finite diff'erence scheme. In this section we apply a finite difference 
scheme proposed by Kushner |23 to the TV-dimensional non-convex equation (|1.1() 
where F is given by H3.1|l and the coefficient a is independent of x. 

We will assume that (Al) and (A2) of SectionOland that the following assump- 
tions hold: 

(A4) a" is independent of x, 

(A5) aS-E,^JaSl>0, i^l,...,N, 

(A6) Ell {< - E,^. + l^f (^)l} < 1 in 

(A7) sup„ {inf^c" - 2VA[6"]i} =: Aq > 0. 

(A8) (i) g e C"'i(R^) or (ii) g G C"^i(M^) and \D^g-\o < C. 

Here we need (A4) in order to prove condition (S4) of Section|21 for more on this 
see To avoid (A4) we must use the much more difficult methods of |3| or PU] . 
We will not consider this here. Condition (A5) simply says that a is diagonally 
dominant. This is a standard condition [23 and implies that the scheme (IB.ip 
below is monotone. Conditions (A6) is a normalization of the coefficients in (|A.lll . 
We can always have this assumption satisfied by multiplying equation IjA.ip by 
an appropriate positive constant. Conditions (A7) and (A8) together with (Al) 
and (A2) assure that the solutions of the various schemes (e.g. (|1.3|l 'l belong to 
C'''^(IR^). Under these assumptions the solutions of various equations (e.g. Hl.l|l l 
will also belong to C'''^(M^). We refer to the Appendix for the proof of these facts. 
Condition (A8) is a regularity condition on cf. (C3) and (S7). 
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The difference operators we use are defined in tlie following way 
A^.w(x) ~ ± — {w{x ± ah) — w{x)}, 
Al.w(x) — —7t{w(x + Cih) ~ 2w(x) + w(x — Cih)}, 
A^.^.w{x) = 2/i2^'^^^^-' ^ ~^ ^^^^ ^ ^ ~ ^^^^^ 

— 2^{2u'(a;) + w(x + eih — ejh) + w{x — ah + ejh)}. 

Let — max{6, 0} and — (—6) + . Note that 6 = 6+ — 6^. For each x, t, pf, 
Au, Af^, i,j = l,...,N, let 

Fix,r,pf,Au,Af^) 

"s-^ j=i 

- bt+ {x)p+ + bf- (xK] + (x, r) } , r - .g(x) I . 

Now we can write the finite difference scheme in the following way, 

(5.1) F{x,Uh{x),A^.Uh{x),Al^Uh{x),A^^^^Uh{x)) = 0. 

This is a consistent and monotone scheme. 

In order to get our result, we must define S in (|1.3|) and prove that conditions 
(A3), (SI) - (S7) of Sections 01 and 0] hold. We have moved most of the details 
to Appendix ^ where a more general problem is considered. To see how S may 
be defined, see (|B.3|I in Appendix IbI Condition (SI) holds by monotonicity of the 
scheme, (S2) holds trivially, (S3) holds with following estimate: 

\F{x,v,Dv,D^v) ~ S{h,x,v{x), H^)| < K{\D^v\oh + {D'^vloh^), 

for any v G C"'(R^). Condition (S5) holds by "convexity" of the sup-part of the 
scheme, (S6) holds trivially, and by (A8) we immediately get (S7). The only difficult 
condition is (S4). To prove it we need very precise a priori estimates on the scheme 
provided by Lemmas IB. II and IB. 21 in Appendix Note in particular that the 
bounds in (S4) are independent of the penalization parameter e. 
In view of Theorem 13 .51 we have the following result: 

Proposition 5.1. Assume (Al),(A2), (A4) - (A8) of Sections\^and\^hold, u is 
the viscosity solution of p. II) and Uh is the solution of (|5.1|l . Then if h > is 
sufficiently small, 

\u-Uh\Q <Ch^/^ in R^. 
If in addition \D'^g~\o < oo, then 

\u-Uh\o <Ch^/'^ in R^. 
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In the convex case under similar assumptions the rate is 1/2 when a is inde- 
pendent of X and at least 1/5 in the general case 0. For one-dimensional 
non-convex problems the rate is at least 1/5 jl4j . and for first order problems the 
rate is again 1/2 9 . 

5.2. Control schemes. In this section, we consider a so-called control schemes 
introduced in the second order case by Menaldi The scheme is defined in the 
following way, 

(5.2) Uhix) = min |(1 - hc^ ix))Utuh{x) + hfix)}, 

where 11^ is the operator defined by 

nUix) ^ 

— J2 {^i^ + hb^'ix) + Vhaiix)) + 0(x + hb'^ix) - VhaUx))) , 

m—l 

and (jf^ is the m-th column of . In the convex case a fully discrete method is 
derived from (|5.2|l and analyzed in [7]. The authors also provide an error bound 
for the convergence of the solution of the fully discrete method to the solution of 
the scheme H5.2|l . 

In this case we only need to assume conditions (Al), (A2), (A7), and (A8), 
in particular a may depend on x. All condition (SI) - (S8) then holds, and the 
consistency condition (S4) takes the form 

\F{x,v,Dv,D^v) - Sih,x,vix), H^)| < Ki\D\\o + \D\\o + \D''v\o)h, 

for any v g C"'(R^). We refer to ^ for the proof of these conditions and the precise 
definition of S. The only difficult point is again (S4). To prove this condition one 
must modify the arguments of in a similar way to what we did in the Appendix 
for the FDM. We omit the details. 

In view of Theorem 13.51 we have the following result: 

Proposition 5.2. Assume (A1),(A2), (Al), and (A8) of Sections\^and\^hold, u 
is the viscosity solution of (|l.l|l and uu is the solution of (|5.2|l . Then if h > is 
sufficiently small, 

\u-Uh\o<Ch^/^^ m R^. 
If in addition |_D^g~|o < oo, then 

\u-Uh\Q <Ch^/^ in R^. 

In the convex case under similar assumptions the rate is at least 1/4 [T^, if the 
solution in addition has 3 bounded derivatives then the rate is at least 1/2 j24j . 
For one-dimensional non-convex problems the rate is at least 1/10 ^1], and for first 
order problems the rate is again 1/2, see e.g. 0. 

Appendix A. Estimates on the Isaacs equation. 

In this section we will give well posedness results and very precise a priori bounds 
for the Isaacs equation 

(A.l) inf sup |-tr \a°'''^ (x)D^u] - b"'^\x)Du + c"^''(x)u - f"''^(x)] = 
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in M^, where a = aa^ for some matrix (function) a. We take the fohowing as- 
sumption: 

(Bl) c > and there is a constant C independent of a, (3 such that 

We start by existence, uniqueness, and L°°-bounds on the solution and its gra- 
dient. 

Lemma A.l. // (Bl) holds and sup„_^ {inf^ c"-'' - [cr"-^]? - [6"^^] i } > 0, then 
there exists a unique solution u of (|A.1|I satisfying the following hounds: 

""""^^jhs:^' '^--'"^iT inf.c-./^-[^"-^]f-[^"-^]i - 

Remark A.l. Usually the assumption on c is c > A > sup^ ^ + 
and all estimates are given in terms of A instead of c. For our purpose this is not 
good enough, since we need to consider limit problems where for some values of 
a,/3, both I/I and |c| blow up, while for others they both remain bounded (cf. the 
penalization method). 

Proof. Existence and uniqueness follows from the (strong) comparison principle and 
Perron's method Let 

M := sup— — g, 

then the first bound on u follows from the comparison principle after checking 
that M {—M) is a supersolution (subsolution) of l|A.l(l . To get the bound on the 
gradient of m, consider 

TO := sup {u{x) — u{y) — L\x — yW . 

If by setting 

i«io[c"''^]i + [r-^ii 



we can conclude that m < 0, then we are done. Assume for simplicity that the 
maximum is attained in (x, y). 11 x = y then m = and we are done. If not, then 
L\x — y\ is smooth at (i, y) and a standard doubling of variables argument leads to 
TO < 0. Since the maximum need not be attained, we must modify the test function 
in the standard way. We skip the details. (The interested reader can have a look at 
the appendix of where the above argument is given for a linear equation.) □ 

Now we proceed to obtain continuous dependence on the coefficients. Let u solve 
the following equation: 

(A.2) inf sup {-tr [a°'''^ {x)D'^u] - b°'^'^[x)Du + c"^'3(2:)u - /"^^(x)} = 
in M^, where a — aa'^ for some matrix (function) a. 

Lemma A.2. If u and u are bounded Lipschitz continuous solutions of (|A.1(I and 

(|A.2|I respectively, and that both sets of coefficients satisfy (Bl). Then 

\u - u\o < sup— _ a - (To 

a^f3 mt^cV mt^c 
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+ sup . I j 2L\b ~ b\o + M\c - c|o + 1/ - /loj, 

where L = [u\i V Af = \u\q V |u|o, a?id 

if2 = 32Lsup|4i[a]2 A [a]? + 2L[6]i A + M[c]i V [c]i + A 

Outline of proof . Define 

m sup|w(x) - u{y) - ~ - + I . 

x.,y \ d J 

Then do doubling of variables using the 3 last terms in the above expression as test- 
function. Using the definition of viscosity solutions and subtracting the resulting 
inequalities lead to 



< sup I - tr[a(y)y] + ti[a{x)X] - b{y)px + b{x)py 

a, 13 ^ 

+ c(y)w(y) - c{x)u{x) - f{y) + fix)y 



where a;,y is the maximum point for m and {px,X), (—py,Y) are the elements 
in second order semi-jets in for u, u given by the maximum principle for semi- 
continuous functions Now we note that by using Lipschitz regularity of the 
solutions, a standard argument yields 

\x-y\<SL. 

So using Ishii's trick \l'6l pp. 33,34] on the 2nd order terms, and a few other 
manipulations, we get 

< snp{-\aix) ~ a(y)|2 + 2L\b{x) - b{y)\ + Ce{\ + \xf + \y\^) 

+ M\c{x) - c{y) I - (inf c V inf c)m +\f{x)~ f{y) \ ] . 

Some easy manipulations now lead to an estimate for m, and using the definition 
of m, we obtain an estimate for \u — u|o depending on 5 and e. We finish the proof 
by minimizing this expression w.r.t. 5 and sending e — ^ 0. □ 

For more details on such manipulations, we refer to |16II17| . 

Appendix B. Estimates on a finite difference scheme 

In this section we apply a finite difference scheme proposed by Kushner 12 Ij to 
the iV-dimensional Isaacs equation ljA.l|) with coefficient a independent of x. 
We will assume that (Bl) of Section fXl and the following assumptions hold: 



(B2) a is independent of x 

'y I 



(B3) oflf ~Y..,^^Wtf\>^. z = l,...,7V, 



in R^. 



(B4) Ef=i {«r/ - Y.,^. Kf\ + \bt''{x)\} < 1 

Let us define the scheme. For each x, r, pf, An, Afj, j, j = 1, . . . , A^, let 
F{x,r,pf,Au,Af.) 
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N a,l3 a,0- 



1=1 J^l 



Using the difference operators A^. , A^., A^.^^. defined in Section [S] we can now 
write the finite difference scheme in the following way, 

(B.f) F{x,uh{x),l^%uh{x),^l^uh{x),^%^^uh{x)) = 0. 

This is a consistent and monotone scheme. In the following it will be convenient to 
use an equivalent formulation of this scheme (see [3] for more details): 

(B.2) Uh{x) = inf sup | ,2aB( \ 

J2 p'^-'^i^, ^ + ^>hix +z) + h^r-f'ix)) \, 



where 



i=l j^i 

p-'^x, X ± e.h) = ^ _ ^ ^ + hbr^H-), 

a.fi± 

a - 

p°''^{x, X + Cih ± Cjh) = — , 

a ■ 

p"''^(x,a; - eih ± ejh) = — , 

and p°'''^{x,y) = for all other y. Let h < 1. Note that by (B3) and (B4), 
< p"-'^{x,y) < 1 for all a,P,x,y. Furthermore J2zehz" P""'^ i^^ ^ + z) = 1 for all 
a, P, X. 

For the readers' convenience we will state explicitly the function 5* (as in 
and H3.2|l ') corresponding to this scheme: We set '■— <j>{x + •) and 

(B.3) 5(/i,2;,r,[0]^):= 



1 



inf sup , , „ 



+ c"'^(x)r-r''5(y) 



E p'^-'iy,y + z)mz)-t 

We use fix point arguments to prove existence, uniqueness, and a priori bounds 
for equation (|B.1|I . 

Lemma B.l. Assume (Bl) - (B4) hold and 

sup (inf c"''^ - 2\/lv 1 1 =: Ao > 0. 

Then there exists a unique solution Uh G C°'^(]R"'^) of the scheme (jB.l|l satisfying 
the following bounds 

I I - / r -I ^ l+h^inf^C.I^ 1^ \l + U Jl 

u/io<sup— — T, u/,.i<sup- 



inf, c"'^ ' ' - inf, c"./3 - 2V^[6"''3]i 
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Proof. Define Th : Ch(R^) ^ Cb(R^) in the following way: 
Thvix) :— inf sup I —5 ^ 

For u,v G Cf,(M^), we subtract the expressions for T/jU and ThV. After we use the 
inequality inf sup(- • • ) — inf sup(- • ■ ) — supsup(- •• — •••), the properties of p"^'^ , 
and (Bl), we obtain 

Thu{x) - Thv{x) 

- ^'^P 1 1 ^ ^2 tf ro^ P 51 P"''^{X,X + Z)\U[X + Z)-V{X + Z)\\ 

Since we may reverse the roles of u and v, we see that Th is a contraction in 
(Cb(]R^), I • |o)- Banach's fixed point theorem then yields the existence and unique- 
ness of a u/i g Cb(R^) solving ljB.21) (and IjB.ip ). The estimate on \uh\o follows 
easily from the identity |u;i|o — \ThUii\o- 

We proceed by proving that Uh has a bounded Lipschitz constant assuming 
for simplicity that c"'^ is independent of x. Let v G C°^^(M^) and subtract the 
expressions for Thv(x) and Thv(y): 

Thv{x) - Thv{y) < 

sup I -I , A aa i \p°'''^{x,x + z)iv{x + z)-v{y + z)) 

+ viy + z){p'^'^ix, x + z)- p^-^{y, y + z))] + h^ir-^ix) - r^iy))^ |. 

In the right-hand side the first sum is bounded by [f]i|a; — and by using the 
definition of p"'*^, the second sum is equivalent to 



N 



i=l 

< 2VNh''\b''-l^{x) - b''-l^{y)\[v]i = 2VNh^[b'''%[v]i\x ~ y\. 

Let C"'^ := 2\/]v'ft,^[6"'^]i. By the above expressions, and by exchanging the roles 
of X and y, we obtain the following estimate 

(B.4) \Thvix) ~ nviy)\ < 

Hence T^v G C°''^{R^), and Uh G C°'i(K^) since Uh = lim^^ooiThYvo for any 
vo G C'^'^(R^). Furthermore, since c"'^ > C"'^ -I- Aq the estimate on follows 
easily from the identity [uh,]i = [ThUh]i and HB.4|I . 
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When c depend also on x we obtain an expression like (|B.4(I with c"'^ and 

snp^Af"^"]! replaced by inf, c"-'' and sup„_^ + [c"-^]i 'it+'^iHf ) re- 

spectively, and hence the lemma would hold again. □ 

Using a standard maximum principle type of argument, we now derive a priori 
estimates on the continuous dependence on the data. 



Lemma B.2. Assume (Bl) - (B4) hold anduh,Uh e C°'_i(R^). If Uh solve (|R1T) 
with data (a,b,c,f) and Uh solve (jB.l|) with data (a,b,c,f) (same a!), then 

\uh ~ Uh\o < sup — 1. _ \2L\b - b\o + M\c - c\o + |/ - /|o|, 

a^p mix c V mt^ c L J 

where L = VN[uh]i V [uh]i, M = \uh\o V \uh\o. 

Proof. We will assume that sup(u — u) — {u — u){x) > 0. The general case follows 
from standard modifications to the proof below. 

Since the scheme (jB.l|l is monotone, at the maximum point x we have 



N a 13 a, 13+ a, 13- 

Y: [ - + E ( - + K - u,){x) < 0, 

i=l j/i 



and 



N 



1=1 

At the point x, we subtract the equations for u/j and Uh. After some rearranging 
using monotonicity of the scheme (the above two inequalities) we get 

< sup E - + (K - &r)(^)A- 1 u^,ix) 

- c{x){uh - Uh){x) - ui,{x){c- c){x) + (/ - f){x) 

This (almost) immediately gives the upper bound on — Uh- Reversing the roles 
of Uh and Uh gives the lower bound and the proof is complete. □ 
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